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The B*B7Z coupling is a fundamental parameter of chiral effective Lagrangian with heavy-light 
mesons and can constrain the B ^ nlv form factor in the soft pion limit which will be useful 
for precise determination of \Vi,h\- We compute the B*Bn coupling with the static heavy quark 
and the 0(a)-improved Wilson light quark. Simulations are carried out with «/ = 2 unquenched 
12^ X 24 lattices at j3 = 1 .80 generated by CP-PACS collaboration. Following the quenched study 
by Negishi et al., we employ the all-to-all propagator with 200 low eigenmodes as well as HYP 
smeared link to improve the statistical accuracy. 
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1. Introduction 

One of the major subjects in particle physics is to determine the CKM matrix elements in order 
to test the standard model and find a clue to the physics beyond. Among all the components, \Vub\ 
is very attractive, since it is determined from the electro-weak (EW) tree-level processes and hence 
gives a good reference point for the standard model prediction. 

Despite its importance, \Vub\ is known only with 10% accuracy from inclusive decay and 20- 
30% accuracy from exclusive decay. Although the inclusive determination has been developed for 
the last decades, it suffers from sizable higher order correction in OPE due to the limited kinematic 
range to avoid the background from B — s- Xclv. In order to reduce the error, significant theoretical 
improvement is needed. On the other hand, the determination from the exclusive B — > Tilv decay 
is more promising, owing to the significant progress in the experiment. Recently BABAR collabo- 
ration gave a very precise determination of the form factor up to \Vub\ [Hjl- The only problem is the 
theoretical uncertainty in the form factor. In view of the present and future experimental situation, 
in order to determine \Vub\ one only needs to know the form factor at any single point in in 
principle. In practice, what would be the best choice of q^l Normally one can choose the for the 
smallest nonzero recoil momenta. However, there is an alternative choice of q^ ~ m\, where one 
can use the symmetry relation. In this choice the B nlv form factor f^{q^) can be expressed in 
terms of the B* meson decay constant fs* and the B*Bn coupling gb as 



where v is the velocity of the B meson, k is the pion momentum and A = mg* — mg. Therefore one 
can reduce problems of computing the form factor to simpler problems of computing the decay 
constant and effective coupling. Since the B* state is below the Bk threshold, it is not possible to 
measure the B* —>■ Bn decay experimentally and only lattice QCD can provide precise information 
of the B*Bn coupling. 

In lattice simulation, the heavy-light decay constant has been studied extensively and is ex- 
pected to be determined precisely with nonperturbative accuracy in unquenched QCD near future. 
However, there has not been much progress in the study of the B*Bn coupling, in particular for 
unquenched QCD. Therefore, in this report we present our recent study of high precision determi- 
nation of the B*Bn coupling in unquenched QCD. 

At present, one of the promising approaches to computing the B*Bn coupling is to use the 
nonperturbative heavy quark effective theory(HQET) including 1 /M corrections. However, it is 
very difficult to calculate the matrix elements for heavy-light systems with HQET. This is because 
in the heavy-light system the self-energy correction gives a significant contribution to the energy so 
that the noise to signal ratio of the heavy-light meson correlators grows exponentially as a function 
of time. In fact, recent results of ^oo are 



these have about 5% and 15% statistical errors for quenched and unquenched cases, respectively. 
But such accuracies will not be sufficient to test new physics. Therefore significant improvements 
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for statistical precision in HQET are needed. Fortunately the two techniques to reduce the statistical 
error are developed recently, which are the HYP smearing |0] and the all-to-all propagators with 
low mode averaging Calculation of goo in quenched QCD using these two methods was carried 
out recently by Negishi et al. It is found that the statistical accuracy is drastically improved as 

|.o = 0.517(16),tat. for ?i/ = 0. (1.4) 

Our ultimate goal is to extend the above strategies to unquenched simulation and give precise 
values of the B*Bti coupling with 2 + 1 flavors in the continuum limit. In this report, we present our 
high precision study of the static B*Bu: coupling in ?ij = 2 unquenched QCD combining two tech- 
niques of the HYP smeared link and the all-to-all propagators. Our results also suggest that these 
techniques can be useful to precisely calculate other physics parameters for heavy-hght systems. 



2. Lattice observables 

The B*Bn coupling can be obtained from the form factor at zero recoil which corresponds to 
the matrix element 

(B* {pB* , A ) \Ai \B{pb)) |^=;i3=o = {^B + )A i {q^ = 0) , (2.1) 

where Ai{q^ = 0) is the matrix element of B to B* at zero recoil with axial current Aj = YJsJiY 
and A stands for polarization [^. In the static limit, 

g^=Ai{q^=0) (2.2) 

holds, so that goo can be evaluated by lattice calculations from the ratio of the 3-point and 2-point 
correlation functions as 

«-«('-'^>-|^. <") 

where C2{t) = (0|^b(0^b(0)|0), €3(1, Ia) = {0\ffB;{t + tA)Ai{tA)ffl{0)\0) are 2-point and 3-point 
functions and ^b, ^^b* are some interpolation operators for the B meson and B* meson with polar- 
ization in the i-th direction. The lattice HQET action in the static limit is defined as 

l + 7o 



S = Y^m—-^ h{x)-Ul{x-4)h{x-4)) 



(2.4) 



where h{x) is the heavy quark field. The static quark propagator is obtained by solving the time 
evolution equation. As is well known, HQET propagator is very noisy, and it becomes increasingly 
serious as the continuum limit is approached. In order to reduce the noise, the Alpha collaboration 
studied HQET action in which the link variables are smeared to suppress the power divergence [Q]. 
They found that the HYP smearing significantly suppresses the noise of the static heavy-hght me- 
son. The error is further suppressed by applying the all-to-all propagator technique developed by 
the TrinLat collaboration [||]. We divide the light quark propagator into two parts: the low mode 
part and the high mode part. The low mode part can be obtained using low eigenmodes of Hermi- 
tian Dirac operator. The high mode part can be obtained by the standard random noise methods 
with time, color, and spin dilutions. Combining these propagators, we can obtain the 2-point func- 
tions for the heavy-light meson which are averaged all over the spacetime. Similarly, the 3-point 
functions can be divided into four parts : low-low, low-high, high-low and high-high parts. 
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3. Simulation details 

Numerical simulations are carried out on 12^ x 24 lattices with two flavors of 0(a) -improved 
Wilson quarks and the Iwasaki gauge action at j3 = 1.80 corresponding to a^^ = 0.9177 (92)GeV. 
We make use of 100 gauge configurations provided from CP-PACS collaboration through JLDG [^]. 
We use the 0(a) -improved Wilson fermion for Ught valence quark with clover coefficient Csw = 
1.60. In the HQET, we use static action with HYP smeared links with smearing parameter val- 
ues {a\,a2,(X'i) = (0.75,0.6,0.3). The B and B* meson operators are smeared with a function 
^{r) = exp(— 0.9r). We obtain the low-lying eigenmodes of the Hermitian Dirac operator using 
implicitly restarted Lanczos algorithm. The low mode correlation functions are computed with Nev 
= 200 eigenvectors of the low-lying eigenmodes of the lattice Hermitian Dirac operator. The high 
mode correlation functions are obtained by using the quark propagator with the source vector gen- 
erated by the complex random noise. The number of the random noise and the number of time 
dilution for each configuration are set to A^;- = 1 and Nf^ = 24, respectively. This setup is based on 
the experience from the work by Negishi et al. [^. 



4. Results 

We show the result of the 2-point correlator with K = 0. 1430 in Fig. |l|, which displays whether 
the improved technique works successfully. Indeed, both the low mode part and the high mode 
part have small statistical errors. The effective mass of the 2-point function is shown in Fig. ^. We 
observe a nice plateau at f > 4. From this result we take f^=5 as a reasonable choice of the time 
difference between the current A,- and the B meson source. 



k=0.1430, smeared ,100configs.Nev=200 




k=0.1430,smeared,100configs,N =200 
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Figure 1: Low mode and high mode contribu- 
tions to the 2-point functions for K = 0.1430. 
Blue and red symbols represent low and high 
mode contributions, respectively. We can see 
low mode becomes dominant for t larger than 
five. 



Figure 2: The effective masses of the 2-point 
and the 3-point functions with K = 0.1430 and 
tA=5. Blue line represents the fit for the 2-point 
functions only. Red line corresponds to simulta- 
neous fit for the 2-point and 3-point functions. 



Fig. H shows the time dependence of the 3-point functions for K = 0.1430. Fluctuations of 
high-modes for 3-point functions are indeed suppressed, so we can get a good plateau for the 
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k=1430 , smeared ,100configs,Nev=200 




Figure 3: Contributions of low-low, high-low, 
low-high and high-high parts to the 3-point func- 
tions for K = 0.1430. 



K=0.1430, smeared, 1 0Oconfigs.N, =24 




Figure 4: The ratio of the 3-point and 2-point 
functions at each t for K = 0.1430. The solid 
line represents the ratio of fit parameters Z3 /Z2. 



effective mass of 3-point functions as in Fig ^. We use the following fit functions 

C2it)=Z2exp{-mt), ^^^^ 
€3(1) =Z3exp(— mf), 

where Z2 and Z3 are constant parameters and m corresponds to the heavy-Ught meson mass. The 
bare B*B7i coupling can be obtained by the ratio of the fit parameters as = Z3 /Z2 . Fit ranges 
for the 2-point and the 3-point functions are 5 ~ 10 and 8 ~ 10, respectively. Fig. ^ shows that the 
effective masses and the fit results are all consistent as we expected. 

Fig. ^compares the results of the B*Bn coupling for K = 0. 1430 determined from two different 
methods. It is found that the ratio C3{t)/C2{t) at each t and the ratio of the fit parameters Z3/Z2 
give consistent value of the B*B7i coupling. We will use the latter result Z3/Z2 to determine gco in 
the following analyses. The physical value of the B*B7i coupling is obtained by multiplying the 
renormalization constant. We use the one-loop result of renormalization factor for the axial vector 
current 

Ai = IkuoZa (l+bA-) Af, (4.2) 

where the gauge coupling = 3.155 and Za = 0.932 for jS = 1.80 as given in Ref. We 

arrive at our preliminary results of for our k values in Table 1. 



K 


0.1409 


0.1430 


0.1445 


0.1464 


8°° 


0.612(5)stat 


0.598(5)s,at 


0.591 (4),tat 


0.578(5)stat 



Table 1: Preliminary results of g, 
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We take the chiral extrapolation of the B*Bn couphng using the data with four quark masses. 
Employing the following fit functions (a), (b), (c), 

(a) g{ml)^ = g{p)+A^ml, 

(b) g{ml)2 = g{0)+Aml+A2{mlf, 

(.) .K),.„o)(,-„o)^^|l„gK))..,.i...,.i)^ 

we carry out the linear extrapolation, the quadratic extrapolation, and the quadratic plus chiral log 
extrapolation where the log coefficient is determined from ChPT[Q]. We use three, four and four 
data point for the fits, respectively. We obtain physical values of the B*Bn coupling in the chiral 
limit as goo = 0.57(1), 0.57(2), 0.52(1) from the linear fit, the quadratic fit and the quadratic plus 
chiral log fit, respectively. We take the average of the results from the linear fit and the quadratic 
plus chiral log fit as our best value and take half the difference as the systematic error from the chiral 
extrapolation. Other systematic errors are the perturbative error of 0{a?), and the discretization 
error of C?((aA)^). Including these errors estimated by order counting, our preliminary result of goo 
is 

= 0.55(l)stat.(3)chiral.(3)pert.(6)disc. at j8 = 1.80. (4.3) 

We find that the discretization error is dominant in our simulation on this coarse lattice. 
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Figure 5: Comparison with other calculations [|[ |[ |^]. In our results, small and large errors represents 
statistics error and perturbative error respectively. 



5. Conclusion 

In this report, we computed the B*Bn coupling on unquenched lattices using the HYP smearing 
and the all-to-all propagators. Using the low mode averaging with 200 eigenmodes, the statistical 
error becomes tiny for all the quark masses, giving ~ 2% in the chiral limit. However, since 
the dominant error is from the discretization for our simulation on this coarse lattice, we need to 
simulate on finer lattices. In Fig. ^, we compare the recent results of the B*Bn coupling [^, ^, |6|]. 
The improvement in statistical precision is drastic, which proves the power of the improvement 
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techniques examined in this report. This result impUes that one can also precisely calculate other 
quantities such as the B meson decay constant or g with 1/M corrections. 
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